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Structural and magnetic properties of the doped terbium manganites (Tb,d)Mn03 {A = Gd, Dy 
and Ho) have been investigated using first-principles calculations and further confirmed by subse¬ 
quent experimental studies. Both computational and experimental studies suggest that compared 
to the parent material, namely, TbMnOa (with a magnetic moment of 9.7 ps for Tb®+) Dy- and 
Ho- ion substituted TbMnOs results in an increase in the magnetic moment lO.Gps for Dy®"*" 
and Ho®"*"). The observed spiral-spin AFM order in TbMnOs is stable with respect to the dopant 
substitutions, which modify the Mn-O-Mn bond angles and lead to stronger the ferromagnetic com¬ 
ponent of the magnetic moment. Given the fact that magnetic ordering in TbMnOs causes the 
ferroelectricity, this is an important step in the field of the magnetically driven ferroelectricity in 
the class of magnetoelectric multiferroics, which traditionally have low magnetic moments due to the 
predominantly antiferromagnetic order. In addition, the present study reveals important insights 
on the phenomenological coupling mechanism in detail, which is essential in order to design new 
materials with enhanced magneto-electric effects at higher temperatures. 


The phenomenon of inducing electric polarization in 
a material by the application of an external magnetic 
field, or vice versa, has attracted considerable attention 
from the point of view of fundamental science challenges 
as well as their envisaged appealing poten tial for appli¬ 
cations in the general area of spintronics.^i^ Since the 
discovery of TbMnOs,^ in which ferroelectric polariza¬ 
tion (~600/rCcm“^) is induced at its magnetic transition 
^27 K and Neel temperature is observed ^42 K, the 
orthorhombically distorted perovskite rare-earth man¬ 
ganites, RMnOa, have become the most notable class 
of materials due to their high magnetoelectric coupling 
(i.e., magnetically induced ferroelectricity).® However, 
due to the relatively low value of electric polarization 
100/iCcm“^) and low Neel temperature 27K), 
these materials are not suitable for most technological 
applications 

Owing to the intrinsic capability of the perovskite 
structure to host ions of different sizes, chemical doping is 
considered as a new knob to tune the physical properties 
of this class of materials.® Both computational and exper¬ 
imental techniques hav e been used to study a variety of 
dopants in perovskiteIn the particular case of 
TbMnOs, dopants with different ionic size and different 
oxidation state such as Ru,^ Ca, b^B0 | SrpH! Na,l^ Al,l^ 
Sn,l23lLa,[^Dy,l2llY,l^ and HcPUm have been introduced 
at the Tb site to tune the magnetic properties. On the 
other hand, dopants with small ionic size such as Cr,l^ 
Sc^ Co,m or FtP have been introduced at the Mn 
site to increase the magnetic transition temperature and 
to induce ferromagnetic (FM) interactions. 

At ambient conditions, the perovskite RMnOs crys- 
talizes in either orthorhombic or hexagonal phases. 
With large ionic size R (such as in the cases of Pr, Nd, 
Tb and Dy) the orthorhombic phase {Phnm) is found 
to be stable (Fig 1(a)), while with small ionic size R 
(such as in Ho, Er, Yb, and Lu) RMnOa stabilizes in the 



FIG. 1: (Golor online) (a) Grystal structure of perovskite 
rare-earth manganites TbMnOs consists of corner sharing 
MnOe octahedra, with the Mn ions in the center of the octa¬ 
hedral. The Tb and dopant A atoms are at the center of a 
cube formed by MnOe octahedra. Green, blue, pink and red 
spheres represent the Tb, A, Mn and O atoms, respectively, 
(b) DFT based energetics comparing two type of magnetic 
ordering considered in (Tb,d)MnOs (where A = Gd, Dy and 
Ho) namely E-AFM type and spiral-spin ordering. Gonsid- 
ered both E-type AFM spiral-spin ordering are also shown. 
For all cases spiral-spin ordering is energetically found favor¬ 
able, hence as a reference is set to be zero. 

hexagonal phase {P63cm, at room temperature). The 
orthorhombic RMnOs with small R ions (Ho and Lu) 
can also be synthesized in a metastable orthorhombic 
phase by high pressure synthesis techniquP^ or special 
chemical methods.l^ Owing to the several magnetic tran¬ 
sitions occurring at low temperatures, RMnOa, (i?=Tb, 
Dy, Lu) with orthorhombic distorted perovskite structure 
have recently been identified as candidates for magnetic 
refrigeration utilizing the magnetocaloric effect (MCE).® 
However, due to the relatively broad temperature depen¬ 
dent peak in magnetic entropy change, the under-lying 
phenomenon governing the origin of such a large mag¬ 
netocaloric effect (MCE) is still unclear. Since, one can 
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treat the magnetic structure and lattice distortions si¬ 
multaneously in first principles computations, these tech¬ 
niques can be used to understand and explain the phe¬ 
nomenological coupling mechanism and can provide nec¬ 
essary insights to design new materials with enhanced 
MCE operateable at higher temperatures. 

In this contribution, using first-principles density- 
functional theory (DFT) based methods we have carried 
out a careful investigation of structural, electronic, and 
magnetic properties of pure and doped terbium mangan- 
ites (Tb,^)Mn03 {A = Gd, Dy and Ho). First, we have 
compared the relative energetics of two types of mag¬ 
netic ordering anticipated in (Tb,H)Mn03, namely, the 
E-AFM type (observed in HoMnOs) and the spiral-spin 
ordering (observed in TbMn03). Next, we have studied 
the effect of dopant A (where A = Gd, Dy and Ho) on the 
electronic and magnetic properties of (Tb,H)Mn03. Fi¬ 
nally, we have compared our computational results with 
experimentally obtained magnetic properties of pure and 
doped TbMn03 bulk ceramics. The present study pro¬ 
vides an improved understanding of both the spin struc¬ 
ture and magnetization of doped systems and suggests 
that the substitution of rare-earth ions with R-ions of 
magnetic moments directly affects the ordering of the 
Tb^'*' moments and indirectly affects the Mn ordering 
(causing an increased ferromagnetic component). 

The DFT calculations were carried out using the pro¬ 
jector augmented plane-wave basis functions as embodied 
in the Vienna ab initio simulation package (VASP).^^®^ 
The exchange-correlation interaction was treated within 
the generalized gradient approximation (GGA) using the 
Perdew-Burke-Ernzerhoff (PBE) functional.!^ It is worth 
mentioni ng tha t for many multiferroic materials such as 
TbMn03,E^ HoMn03,® and GaMnyOia,!^ the first- 
principles results are very sensitive to the choice of the 
on-site Goulomb energy f/psHUl Qm. tests with different 
U values confirm that for TbMn03 U = 2 eV for the 3d 
states of Mn reproduces the experimental band gap value 
of 0.5 eV.ESl Hence, all the calculations are performed us¬ 
ing U = 2 eV.!^ A Monkhorst-Pack fc-point meshP^ of 
3x3x2 is employed to produce converged results within 
0.1 meV per formula unit. In all calculations spin-orbital 
coupling (SOC) was included. In doped TbMn03 stud¬ 
ies, one Tb atom was substituted with dopants Gd, Dy, or 
Ho, thereby resulting in a dopant concentration of 25%. 

Powder samples of pure TbMn03 and doped 
Tbo.67Ao,33Mn03 {A = Gd, Dy and Ho) were prepared 
by a solution based citrate route. High purity nitrate 
precursors of Tb, Ho, Dy, Gd, and Cr as well as citric 
acid were dissolved in water and mixed in the stoichio¬ 
metric ratios. After drying, the powder samples were 
annealed at 900°C for 2 hours in an O2 environment. 
Structural characterization by X-ray diffraction (XRD) 
was performed using a Bruker D8 diffractometer using 
Cu- Ka radiation. Subsequent Rietveld refinement of 
the diffraction data was accomplished using the FullProf 
Suite package.!^ The background was fitted by a linear 
interpolation between a set of background points with 


refinable height while refining the structural parameters, 
and then removed to obtain an accurate goodness of fit. 
DG magnetization measurements were performed with 
the Vibrating Sample Magnetometer attachment to the 
Evercool Physical Property Measurement System from 
Quantum Design. 

TbMn03 stabilizes in an orthorhombic perovskite 
structure {Pbnm symmetry) consisting of the corner 
sharing MnOe octah edra, w ith the Mn ions in the cen¬ 
ter of the octahedral.l2S*2Z138i (Tb,A) ions are at the 

center of a cube formed by MnOg octahedra as shown in 
Fig. 1(a). This structure is highly distorted because of 
the mismatch of (Tb,A)-0 and Mn-0 bond distances and 
the Jahn-Teller (JT) distortion of the Mn-0 bonds. It 
has been noticed that multiple factors contribute to the 
spin-orientation of a system. Specially, in orthorhom- 
bically distorted RMO 3 {R = rare earth, M = transi¬ 
tion metal) spin-orientation is dominated by the M-ion 
array and the interaction between the R and M sub¬ 
lattices, such as antisymmetric exchange interactions,!^ 
strong couplin^^and anisotropic interactions.!^ Here, as 
mentioned above two types of the possible magnetic or¬ 
dering arrangements (below ferroelectric transition tem¬ 
perature) in (Tb,A)Mn03 are investigated. Our first- 
principles calculation based energetics confirm that both 
pure and the doped terbium manganites favor the spiral- 
spin type magnetic ordering, shown in Fig.l (b). 

To get insights on the interplay between orbital- and 
magnetic- ordering in these manganites we have analyzed 
the density of states (DOS). DFT-I-U-I-SOG computed to¬ 
tal and atom projected DOS of pure and doped terbium 
manganites are plotted in Fig. 2 (a-d). It is evident from 
Fig. 2 that the large portion of the DOS (-6 eV to 5 eV) 
is contributed mainly due to Mn atoms accompanied by 
Tb atoms. First, we focus on the pure TbMn03 (Fig. 
2a), which is found to be an insulator with energy gap 
of about 0.5 eV.!^ The DOS in the conduction band can 
be explained in terms of optical conductivity spectra and 
first peak near ^1.7 eV is mainly an attribute of the first 
optical transition as observed in earlier optical conduc¬ 
tivity spectra measurements.!^ The noteworthy feature 
of the atom projected DOS is that the peak correspond¬ 
ing to the lowest optical transitions shows a slight shift 
toward the Fermi level with increasing dopant ionic ra¬ 
dius (see Fig. 2 (a-d)). 

Next, we examine the partial density of states (PDOS). 
The orbital-projected DOS in Fig. 3 (a-d) also confirm 
that in the vicinity of the Fermi-level, the majority of 
the DOS in the valence band arises from the d-states of 
the Mn atoms, while in the conduction band, the DOS 
have contributions from both Mn and O atoms. The no¬ 
table point is that the highest occupied level shows O 
2p character, while the lowest unoccupied level has Mn 
3d character and mainly dominated by occupation of the 
d(x'^-z^) orbitals (Fig. 3). Specifically in the region just 
above the Fermi level (between 0-2 eV) the DOS arise 
due to the unoccupied Mn (d) states and above 

2 eV the DOS spectrum is dominated by all Mn (3d) 
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TABLE I: Summary of the measured structural and magnetic properties of TbMnOs and Tbo.eTAo.aaMnOa {A = Gd, Dy and 
Ho). Quantities computed using DFT calculations using the GGA+U methods are provided in parenthesis. 


Atomic 

TbMnO.':^ 

Tbn RyGdn ^^^MnO^ 

Tbn 67Dyn..-3.-3MnO,'^ 

Tbn,fi7Hon.r?riMnO,', 

coordinates 

X y z 

X y z 

X y z 

X y z 

Tb/i? {4c) 

0.98 (0.98) 0.08 (0.08) 0.25 (0.25) 

0,98 (0,98) 0.08 (0,08) 0.25 (0.25) 

0,98 (0.98) 0.08 (0.08) 0,25 (0.25) 

0,99 (0,98) 0.08 (0.08) 0.25 (0.25) 

01 (4c) 

0.08 (0.11) 0.50 (0.46) 0.25 (0.25) 

0,08 (0,12) 0.50 (0,46) 0.25 (0.25) 

0.09 (0.12) 0.51 (0.46) 0,25 (0.25) 

0,08 (0.12) 0.50 (0.46) 0.25 (0.25) 

02(4c) 

0.69 (0.70) 0.33 (0.32) 0.05 (0.05) 

0,69 (0.70) 0.33 (0,33) 0,05 (0.05) 

0.70 (0.70) 0.33 (0.33) 0.05 (0.05) 

0,70 (0.70) 0.32 (0.33) 0.05 (0.05) 

Mn (4c) 

0.00 (0.00) 0.00 (0.00) 0.50 (0.50) 

0,00 (0.00) 0.00 (0,00) 0,50 (0.50) 

0.00 (0.00) 0.00 (0.00) 0.50 (0.50) 

0,00 (0.00) 0.00 (0.00) 0.50 (0.50) 

Structure parameters 

a (A) 

5.293 (5.335) 

5.310 (5,341) 

5.302 (5.327) 

5.271 (5.324) 

b (A) 

5.839 (5.884) 

5.844 (5,914) 

5.833 (5.895) 

5.811 (5.892) 

C (A) 

7.402 (7.481) 

7.424 (7,471) 

7.418 (7.453) 

7.378 (7.449) 

‘^l(Mn-02) 

2.168 (2.246) 

2.168 (2,246) 

2.182 (2.243) 

2.127 (2.235) 

^^(Mn-02) 

1,986 (1,981) 

1.986 (1,981) 

1.944 (1.974) 

1.988(1.975) 

^(Mn-Ol) 

1,907 (1,936) 

1.907 (1,936) 

1.920 (1.930) 

1.896 (1.928) 

Mn - 02 - Mn (deg) 

144.13 (144.16) 

143.8 (144.16) 

145.5 (144.25) 

144.8 (144.15) 

Mn ' 01 - Mn (deg) 

153.5 

153.5 

150,0 

153.3 


2,53 

2,53 

1,96 

3,48 

Magnetic Moment ) 

11,31 ± 0.01(11.29) 

9,6± 0.01 (9.29) 

11.29± 0.01 (11) 

11.1 ± 0.01 (10.7) 



-6 -4 -2 0 2 4 -6 -4 -2 0 2 4 


FIG. 2: (Color online) Calculated total and atom projected 
density of states (DOS) of (a) TbMnOa (b) (Tb,Gd)Mn03 
(c) (Tb,Dy)Mn03 (d) (Tb,Ho)Mn03 atoms. In each case, 
both up and down spin contributions to the DOS are shown. 



FIG. 3: (Color online) Calculated orbital projected density of 
states of O and Mn atoms in (a) TbMnOa (b) (Tb,Gd)Mn03 
(c) (Tb,Dy)Mn03 (d) (Tb,Ho)Mn03 . 


orbitals. On the other hand, in the valence band below 
-2 eV joint contribution from O {2p) and all Mn {3d) 
orbitals has been observed. It is also evident from Fig. 3 
(a-d) that the O {2p) states and Mn (3d) states further 
enhance the strong hybridization between the orbital and 
spin order resulting in the magnetic and structural mod¬ 
ulations, consistent with the Jahn-Teller mechanism.!^ 

The XRD 9-29 data for the pure TbMnOa and doped 
Tbo.er^o.ssMnOa {A = Gd, Dy and Ho) samples is shown 
in Fig. 4(a-d). Peaks corresponding to only the or¬ 
thorhombic symmetry (space group Pbnm) are observed 
indicating that all the samples are phase-pure. Rietveld 
rehnement was used to refine the XRD data of all sam¬ 
ples. Measured and DFT computed structural parame¬ 
ters (lattice constants, atomic positions, etc.) are listed 
and compared in Table 1. The lattice parameters for pure 
TbMnOa are found to be in agreem ent with available 
computational and experimental data.l2S12I1281 slight 
change in lattice parameters due to the doping can be 
understood in terms of the ionic size mismatch between 
host (Tb) and dopants (Gd, Dy and Ho). For instance, 
the lattice constant was expected to decrease when Ho^'*' 
of smaller ionic radius (107.2 pm in 9-coordination) is 
substituted in the site of Tb^"*" with larger ionic radius 
(109.5 pm in 9-coordination). For doped terbium man- 
ganites, the average ionic radius of the A-site ion < r^> 
is 108.741 pm and the ionic-size mismatch represented 
by the variance cr^= Ey^ rf - <ryi>^ is 1.17 pm^, where 
Vi corresponds to the radii of the various A-site cations 
and Yi to their fractional occupancies. The lowered aver¬ 
age radius of the rare-earth ion also explains the reduced 
value of the Mn-O-Mn bond angle (in-plane) as listed in 
Table 1, which as previously mentioned has been shown 
to be crucial in determini^ the magnetic properties of 
the rare-earth manganites.l^ 
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FIG. 4: (Color online) The simulated and measured 
XRD 9-26 scans of orthorhombic (a) TbMnOa, (b) 
Tbo.srGdo.ssMnOs, (c) Tbo.srDyo.asMnOa, and (d) 
Tbo.srHoo.ssMnOa. 



Temperature (K) 

FIG. 5: (Color online) (a) The temperature dependent 
zero-field cooled (ZFC) and field cooled (FC) dc magneti¬ 
zation data for (a) TbMnOa, (b) Tbo.erGdo.aaMnOa, (c) 
Tbo.srDyo.aaMnOa, and (d) Tbo.erHoo.aaMnOa bulk samples. 


To validate our predicted results we measure the mag¬ 
netic properties of these powder samples. The temper¬ 
ature dependent magnetization of the present TbMnOa 
and Tbo.ey^o.asMnOa {A = Gd, Dy and Ho) samples 
is shown in Figure 5 in both zero-field-cooled (ZFC) and 
field-cooled (FC) conditions (with applied field of 50 Oe). 
The peak at ~ 7 K in each of the present samples is due 
to the rare-earth magnetic ordering. No such features 
are observed for the paramagnetic to spiral antiferromag¬ 
netic and incommensurate to commensurate Mn mag¬ 
netic transitions for pure and doped terbium manganites 
samples.l^ A bifurcation of the ZFC and FC curves is 
observed at 60 K and 63 K for Tbo.erDyo.aaMnOa and 
Tbo. 67 Hoo. 33 Mn 03 , respectively, suggesting a weak fer¬ 


romagnetic component of the Mn ordering. The high 
temperature (paramagnetic) susceptibility data was fit 
to a Curie-Weiss law: x=C/(T-6*), where x is the sus¬ 
ceptibility, C is the Curie constant, and 9 is the Weiss 
temperature. The effective magnetization are shown in 
Table 1. The negative Weiss temperatures observed in 
all of the present samples indicate antiferromagnetic or¬ 
dering. The effective magnetic moments were calculated 
from the Curie constants. Gd substitution significantly 
reduces the effective magnetization, whereas Dy and Ho 
substitution only slightly reduce the effective magnetiza¬ 
tion as compared to pure TbMn 03 . This is likely due 
to the high magnetic moments of Tb, Dy , and Ho com¬ 
pared to the moderate moment of Gd.^^^ The measured 
magnetic moment is found to be in good agreement with 
our DFT calculated values (see Table 1). 

Above the Tb ordering temperature at 7 K, the Tb 
sub-lattice is expected to exhibit ordering as a conse¬ 
quence of ^Mn-Tb exchange coupling.lS^ Even below 7 K, 
the ordering of Tb is demonstrated to be related to the 
wave vector of Mn ordering in this temperature regime 
as a result of relatively strong Mn exchange field .1^ Tb 
is thought to be a special case of intermediate coupling 
that leads to complex spin ordering behavior with respect 
to Mn sublattice, whereas Dy-Mn coupling is weak and 
Ho-Mn is very strong which leads to the wave vectors 
completely decoupled or coupled (respectively) down to 
zero temperature. Thus, we can expect that the substi¬ 
tution of 33% Tb with another rare-earth ion will signif¬ 
icantly impact the exchange coupling and consequently 
the ordering of the rare-earth moments at low temper¬ 
atures. The modulation of the Mn-O-Mn bond angles 
and ^Mn-A in the present samples may explain a signif¬ 
icant difference in susceptibility and coercivity observed 
in the magnetic field dependent magnetization behavior 
(not shown), with stronger coupling leading to increased 
susceptibility at low applied fields. Gd substitution in 
TbMnOs slightly reduced the Mn-O-Mn bond angle and 
lead to slightly weakened magnetization values. Dy and 
Ho substitution both resulted in higher Mn-O-Mn bond 
angle than that in pure TbMnOs and enhanced ferro¬ 
magnetic component of the magnetization value. Addi¬ 
tionally, the strong Ho-Mn interaction could account for 
the much larger increase in the magnetization value in 
case of Ho substitution, as compared to the weak Dy-Mn 
interaction present in the Dy substituted sample. How¬ 
ever, any modification to the magnetic ordering of the 
Mn moments (he. the weak ferromagnetism observed 
in Dy- and Ho-doped TbMnOs samples at higher tem¬ 
peratures) is expected to be a result of the theoretically 
and experimentally determined change in the Mn-O-Mn 
bond angles that will directly affect super-exchange in¬ 
teractions. 

In summary, we extend this dopant-mediated search 
for new possible multiferroic candidates by introducing 
rare-earth dopants in the parent material TbMn 03 . Us¬ 
ing DFT based computations we predict that at low tem¬ 
peratures magnetic ground state of both pure and doped 
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material systems is spiral-spin type of magnetic ordering. 
The predicted results are in favorable agreement with ex¬ 
perimental measurements on the same material systems. 
Our theoretical and experimental research efforts have 
led to an improved understanding of both the spin struc¬ 
ture and magnetization of doped systems. The major 
findings of this work are that the spiral AFM order ob¬ 
served in TbMnOa is stable with respect to the present 
A-site substitutions. These substitutions do, however, 
modify the Mn-O-Mn bond angles. The higher Mn-0- 
Mn bond angles, caused by the Dy and Ho substitutions, 
correspond with a stronger ferromagnetic component of 
the magnetic moment. This is a crucial finding that the 
magnetic order of TbMnOa, which is known to cause the 
ferroelectricity, can be preserved while the net magnetic 


moment enhanced. This is an important step in the field 
of the magnetically driven ferroelectric class of magne¬ 
toelectric multiferroics, which have traditionally had low 
magnetic moments due to the predominantly antiferro¬ 
magnetic order.l^ In a broad sense, our research has led 
to an improved understanding of both dopant-medieated 
spin structure and magnetization in TbMnOa. However, 
further studies are needed to determine how these present 
substitutions could affect the ferroelectric polarization. 
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